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We study evolution equations and stationary homogeneous solutions for electric and magnetic field 
amplitudes in a ring cavity with fiat mirrors. The cavity is filled with a positive or negative refraction 
index material with third order Kerr-like electric nonlinearities and also magnetic nonlinearities, 
which can be relevant in metamaterials. We consider the degree of freedom of polarization in the 
incident beam. It is found that considering a magnetic nonlinearity increases the variety of possible 
qualitatively different solutions. A classification of solutions is proposed in terms of the number 
of bifurcations. The analysis can be useful for the implementation of optical switching or memory 
storage using ring cavities with non linear materials. 
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I. INTRODUCTION 

During the last decade, progress in the development of composite materials allowed the experimental observation 
of new optical properties, such as negative refraction index [l|, Q. These materials have both negative dielectric 
permittivity and negative magnetic permeability, a property not found in natural materials. Materials with negative 
refraction index (NRM) have many interesting new features, and a great number of applications have been proposed 
based on their novel properties. 

On the other hand, ring cavities filled with nonlinear media have been proposed for information storage and also as 
the optical analog of electrical transistors Q . Conditions for modifying and tuning optical bistability in these cavities 
are being researched [3, @ ■ It is known that electromagnetic wave propagation in a Kerr-type nonlinear material with 
a positive refraction index (PRM) can be described by an order parameter equation of the nonlinear Schrodinger type 
and the same type of equation can be extended to NRM [Ml- 

A recent work 0 shows that arguments previously used to neglect the magnetic response in a PRM do not hold 
for NRM, where both the polarization and the magnetization may take non-negligible values. Even, it has been shown 
[13, [I3| that a composite metamaterial with negative refraction index could develop a nonlinear macroscopic magnetic 
response. This means that although the medium has a negligible magnetic nonlinearity, periodic inclusions produce 
a nonlinear effective response when the wavelength is much larger than the periodicity of the inclusions. 

In this paper we analyze the equations describing the electric and magnetic fields in a rin g ca vity with flat mirrors 
that contains a NRM or a PRM with electric and magnetic nonlinearities. In a previous work [ij, we limited ourselves 
to a linearly polarized incident beam, now we focus on analyzing what happens when the light polarization is taken 
into account, making it necessary the double of equations to describe the system. The main objective of this paper 
is to explore the various types of stationary homogeneous solutions that the system can display, depending on the 
parameters that describe it. We also perform a brief analysis of the behavior of linearly unstable solutions. 

The paper is organized as follows. In section |n] we present the evolution equations for the amplitude of the electric 
and magnetic fields and the reduction to two coupled Lugiato Lefever equations for a ring cavity. In section IIIII we 
analyze homogeneous stationary solutions. The type and number of bifurcations that can occur, depending on the 
parameters, are also analyzed. Finally, in section El numerical integration is performed in order to analyze the 
temporal evolution of the unstable solutions. In section El we present our conclusions. 


II. CAVITY EQUATIONS 

The equations often used for describing the behavior of electric and magnetic fields, at the plane perpendicular to 
propagation of light in a cavity, are of the type of Lugiato-Lefever equation (LL, see [l5j) with generalized parameters. 
LL equation is a simple mean field model, and has been useful for the analysis of pattern formation in cavities 
containing a Kerr medium that is driven by a coherent plane wave [110. 

First, we analyze the problem of free propagation (without mirrors) of an electromagnetic wave in a nonlinear 
material, either PRM or NRM, and then we use the resulting equations to obtain the behavior in the cavity. 
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A. Evolution Equations in the material 


We consider a plane wave of arbitrary polarization and frequency uiq. We assume that the electric and magnetic 
fields are in the plane x-y and the wave propagates in the z axis. The starting point is Maxwell’s equations and 
constitutive relations for the electric displacement, D = eqE + P, and magnetic induction, B = /xpH + /.toM. An 
isotropic metamaterial is analyzed with the third-order Kerr nonlinear response. The nonlinear relationship between 
the polarization of the material and the electric field is Q: 
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where i,j, k and I correspond to Cartesian axes, and take values x, y or z; repeated indexes involve addition. is 

( 3 ) 

a rank 2 tensor describing linear electric behavior in the material and Xe ijki accounts for third order nonlinearities. 

Zharov et al [l3l| analyzed properties of a microstructured metamaterial made of metallic wires and split ring 
resonators (SRR) embedded into a Kerr permittivity material. They found out that, in such a material, the electric 
fields around the SRR could be very strong. In that case, the electric nonlinearity produces an effective magnetic 
nonlinearity fJ-eff- Expanding yeff in powers of the incident magnetic field, to the lowest non constant order, a 
Kerr-like magnetization can be found. Thus, a similar relationship to m is proposed for magnetization and magnetic 
field (with similar definitions for Xm Xm )• 
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Many proposed metamaterials consist of the repetition of a cubic unit lattice which has the same form when it is 
looked from any of its faces (for example, materials usii^ spherical inclusions [l^, spherical inclusions and wires 
M but also cells specifically designed to be symmetric [l^), making optical properties to be the same in any two 
orthogonal directions. If the wavelength of incident light is much larger than the cell length, the material is expected 
to be isotropic. So, we assume that the material is cubic centrosymmetric and isotropic. As the components of the 

(S') 

tensors Xe/m part only on field convolutions, the tensors are defined, without loss of generality, to be symmetrical 
for a large number of indexes and argument permutations, similar to what has been done in Q, Ch. 2, Sect. 2d. 

We apply the classical multiple scales perturbation technique, and we assume that light is quasi-monochromatic 
and can be represented by a plane wave, with frequency wq and wavenumber fcg. It travels thought the z axis, and is 
modulated by a slowly varying amplitude. The amplitude depends on position R = {X, Y, Z) and time T, variables 
that have characteristic scales much larger than the scales given by 1/fco and I/wq. In summary, the fields E and H 
can written as 


E = -hc.c. 

H = -H(R,T)e*('='>^“‘^'>‘)-hc.c. (3) 

where £ and T-L are slowly varying amplitudes. Similar relationships are defined for P and M; the corresponding 
amplitudes can be written as (for the polarization of material see, for example, EH) 

r = 3xL'^ [AE{\£.f + \£yf)£ + BE/2{£^+£^,)£] 

M = 3x^^^[AMm.\^ + \ny\^)'H + BM/2{nl + nl)-H], (4) 

where the bar stands for complex conjugate, Xe/m ~ Xe/m xxxx parameters Ae/m and Be/m are given by 
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FIG. 1: Scheme of the ring cavity. 
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In the previous expressions we are using the notation X^/m abed ('''without explicit dependence) to identify the Fourier 

/ q'\ 

transform of ^ 2 , evaluated in (wo,a;o, -wq). 

We define the wavenumber fc(w) = a;n(a;)/c; k' and k" are the derivatives of fc(w) at w = wq; the refraction index is 
n{uj) = Cr{to) (which takes a negative value when both and are negative [13)) where triio) = 1+Xb^(w) 

and = 1 + Xm (w) are the relative permittivity and relative permeability, respectively. 

The details of the implementation of multiple scales technique in our case, are essentially the same as those described 
in reference Q, Sect. 2k. After this process, we obtain nonlinear Schrodinger equations for the envelopes of electric 
and magnetic fields defined in the plane perpendicular to the z axis. Applying a change of variables to use circularly 
polarized components, \/2£± = ± i£y, \/2'H± = T-Ly =F iHx, and a coordinate transformation given by ^ = z, 

T = t — k'z, the following equations of evolution in the material can be found 
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where the approximation k^ ~ uiQHoHrZoZr, valid for low dispersion fields, is used. For simplicity of notation we write 
/ij.(wo) = Hr and er{uJo) = £r- 

Linear and nonlinear dissipation is neglected, so that Xb/m ^'^d X^e}m valued quantities. This is an often 

used approach in conventional optics, but dissipation could play an important role in metamaterials. Nevertheless, as 
we comment in the next section, there are ways to reduce dissipation in metamaterials even to negligible levels. On 
the other hand, a generalization of the model for small dissipation is also possible H- 


B. Cavity effects 

We consider a ring cavity with plane mirrors and with a nonlinear material as shown in Figure [T] At the input 
mirror, the incident fields are Ain and iJjn, and the reflection and transmission coefficient are and ti respectively; 
Tr is the round trip time and L is the nonlinear material length. The detuning between the incident field and the 
cavity mode is 9 = cj) mod where (j) is the phase accumulated by the wave around the cavity, including the effects 
of reflections in the mirrors. We study the cavity near resonance, so the mismatch is small and e*® ~ 1 + i9. 
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From Maxwell equations, assuming that nonlinearities are small, it is obtained that, in a given medium (different 
from vacuum), the relationship between electric and magnetic fields is approximately 'H± = rj£± while in vacuum the 
relation is 'H± = r]Q£±, where ry and 770 are the inverses of the material and vacuum impedances, respectively. If we 
consider media with similar impedances (ry ~ ?yo), the wave is mainly transmitted, generating a negligible reflected 
wave . The transmitivity coefficients at both ends of the material are t± = 1± , where t+ holds for £ at point 

b (see Fig. [T]) and for H at point a, and t- holds for % ai h and for £ ai a. 

We can obtain an equation of evolution for electric and magnetic fields in the cavity and it can be shown that under 
the above conditions, the proportionality between them is maintained, so the system is well described knowing only 
the electric field. The procedure is analogous to that performed in [l^ . It is important to mention that from now on 
we are looking at cases where transverse spatial dependence can be neglected, i.e. V\£± — 0. We get 
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In order to reduce the number of parameters, an equation with adimensional quantities can be obtained using the 
following change of variables 
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where p = rA^t- 1 so that 1 — p is a small number of the order of P, and £t is a characteristic field given by 


It is useful to define the parameter B as 
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that can be understood as an average of the electric and magnetic constants Be and Bm- In the general case, using 
symmetry arguments, it can be shown that Ae + BeI‘^ = Am + BmI‘2- = 1. In Ref. [1^ (see p. 197), three possible 
cases are mentioned for only electric nonlinearities (i.e. Am = Bm =0): = 0 (B = 0), for electrostriction, that 

does not act at optical frequencies; Be = QAe {B = 3/2) for molecular orientation effects; and Be = Ae {B = 2/3) in 
materials with electronic response far from resonance frequency. In the rest of our work, we will find analytical results 
for any value of B, but we will limit our numerical results to the case 0 < B < 2, which corresponds to situations in 
which Ae/m ^ 0, and where x^e^ li^is the same sign as X^e^■ The motivation of the last choice is a result of Ref. [l3l| . 
where the authors derive an effective Kerr-like nonlinear magnetization for a split ring resonator; it can be written as 
PeS = P-effo + KH^, where K is a complex number related to geometrical factors, which tends to a real quantity with 
the same sign as the electric nonlinearity, when conductivity tends to 00 . Eq. (jS)) finally reads: 


dA± 

dt' 


— Ain±: ~ (1 + z0)A± -1- iq;(|A±|^( 1 — B/2) -|- |Aip|^(l -|- B/2))A±, 


( 12 ) 


where a is the sign of -I- Xs^t-- Let us note that dispersion effects, proportional to fc", are not included, 

neither in Eq. ([8]) nor in Eq. m- The reason is that dispersion generates a term in (1121) proportional to the small 
parameter 1 — p used in the scaling ® and (US, so that it can be neglected with respect to the other terms. 

Considering transverse spatial effects would produce the inclusion of a Laplacian in Eqs. (HI, transforming them 
into Lugiato Lefever equations. Complex conjugate fields describe the same physical system and the corresponding 
equation is equal to (O with 0 and a with opposite signs, therefore the relevant parameter is a0. Considering a 
metamaterial in this analysis justifies the inclusion of the magnetic nonlinearity, that allows a wider range of possible 
values of B. The sign of the refraction index does not appear explicitly in (1121) . but it would appear in a diffraction 
term if we would take into account transverse spatial effects 1^ . 

Up to now, losses have not been taken into account. In the Gigahertz range, it is possible to build an NRM 
with small (and even negligible) imaginary parts of and as was shown in, for example, [^. but for higher 
frequencies the former is not true. Many proposals for reducing losses at high frequencies are being researched, such 
as [261 - 1^ . If dissipation is not strong, Eq. (TT^ can be generalized in a similar way to what has been done in a 
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previous work [l^. Taking fc(u;o) = fco + iki, if kj <C fcg, the term —ki£± has to be added to the right hand side of 
equation d®]), —ki'H± has to be added to the right hand side of equation © and —Lki£± has to be added to the right 
hand side of equation ([8]). An adequate change of variables is easily found replacing (1 — p) by (1 — p + pLkj) in ([9]) 
and (uni, so that m is still valid. 


III. HOMOGENEOUS STATIONARY SOLUTIONS 

Looking for stationary homogeneous solutions in (I12L we obtain 

= [1 + i0 - ia(l - B/2)\A±\^ - iail + B/2)\A^\^]A±. (13) 

From this equation, using that = 1 we can write 

A/„ = (1 + [a© - (1 - 5/2)/+- (1 + B/2)/_]2)/+ 

(l-A)/„ = (1 + [a© - (1 - 5/2) J_ - (1 + 5/2) J+]2)/_ (14) 

where I± = |A±p, and A = is a quantity that measures the polarization of the 

incident beam, related to the ellipticity y by A = cos^(x/2). 

Given a set of parameters {aQ,B,X), solution sets can be obtained. By definition, 0 < A < 1, but, 

for each solution with A > 1/2 there is another solution that can be found exchanging with so, we work with 
1/2 < A < 1. 

It was found that the type of possible solutions varies according to the parameters of the system. We present a list 
of them for linear, circular and elliptic polarization. 


A. Linear polarization, A = 1/2 

An analysis of Eq. (jl4l) indicates that the symmetric solution (/+ = /_) always exists, but there may be more 
solutions. We found, in terms of bifurcations of intensity, a total of 8 possible cases. These 8 possible cases define, in 
the aQ-B plane, 8 regions, that are limited by 4 curves, as shown in Figure [21(a). 

The symmetric solution (regardless if asymmetric solutions exist or not) presents bistability if a© > indepen¬ 
dently of B. The bistability region (i.e. the range of values of Bn for which there are two stable solutions of 1+ ) is 
limited by a pair of saddle-node bifurcations. It has been extensively studied by other authors (SiBBElia)- 

There might be asymmetric solutions /+ > I-, and analogous solutions with < /_. From Eq. a, looking for 
an expression for which is real and positive, depending on J-|_ for the case /_ ^ I+. it is found that; for B > 1, 
there is a pitchfork bifurcation in I*^ where two asymmetric solutions emerge and exist for all Bn > I*n (solutions 2, 
5 and 8 in Fig. El); for B < 1 and a© > 2\/l — B/B, asymmetric solutions appear between two values of Bm in which 
pitchfork bifurcations occur (solutions 3, 6 and 9 in Fig. |31); and for i? < 1 and a© < 2\/l — B/B, only symmetric 
solutions exist (solutions 1 and 4 in Fig. |3l). 

Finally, it may happen that the asymmetric solution presents bistability, i.e., the asymmetric branches (generated 
by a pitchfork bifurcation) have a curling due to a couple of saddle node bifurcations. To find them, it is necessary to 
find the values of /+ on the asymmetric branches where (/+,/_(/+))/9/+ = 0. For example, in Fig. [31 solution 
8, there are 4 values of /+ where this condition is met, and in solution 5 there are no values that verify the condition. 
The limit curve in phase space aQ-B can be analytically obtained. It is found that these bifurcations arise for values 
of a© greater than 2'/ijB. 

It is not possible to have a symmetric solution with two pairs of saddle-node bifurcations. For this to happen, we 
should have 4 points where (/+,/_ = I+)/dI+ = 0. This can not be achieved because Bn is a polynomial of 
degree 3 in when /_=/+, and this is the reason why there is an empty square in Fig. [31 

The possible cases can be organized in a simple way: for a given set of parameters, we define Up as the number of 
values of Bn where there are pitchfork bifurcations and as the number of regions where bistability arises due to 
a couple of saddle-node bifurcations. There may be only the symmetric solution [up = 0), symmetric solution and 
asymmetric solution between two values of incident field (up = 2) or asymmetric solution starting from an incident 
field value (rip = 1). Independently of the number of pitchfork bifurcations Up, there may be rig =0, 1 or 2. 

Examples of typical solutions are shown in Fig. [31 which also indicates whether the various branches are stable or 
not under homogeneous perturbations. In order to classify the solutions with a unique index we assign a number, 
given by A' = 1 -I- rip -I- 3ns, to each possible case. 
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FIG. 2: Types of solutions in phase space 0:0 — B. For A = 1/2 (a) and A = 1 (d), the number in each region is given by 
= 1 + rip + Sris, where rip is the number of pitchfork bifurcations and ria is the number of bistability regions given by a pair 
of saddle node bifurcations. The analytic expressions for the limits between regions are shown in some cases. For elliptical 
polarization (b) and (c), we show a numerical exploration of possible solutions, and we use the number N = 1 + rip + Sris to 
identify regions, where n/ is the number of saddle node bifurcations that correspond to pitchfork bifurcations for the linear 
polarization case. The different filling patterns correspond to values of rip or n'p equal to 0 (diagonal lines), 1 (horizontal lines) 
and 2 (white). 


B. Circular Polarization, A = 1 

From Eq. da, it can be shown that if A = 1, /_ is equal to 0. For a0 < there is only one possible solution, 
and, for a© > there is the well known bistability generated by a couple of saddle-node bifurcations. Possible 
solution regions are shown in Fig. [21(d). These two only possible solutions for circular polarization are equivalent to 
solutions 1 and 4 for linear polarization (see Fig. o. 


C. Elliptical Polarization, 1/2 < A < 1 

As we have seen in the previous subsections, there are many possible kinds of solutions for A = 1/2, but only two 
for A = 1. For 1/2 < A < 1, up to 11 kinds of solutions were found, but the boundaries of them are related to the 
roots of degree 9 polynomials, difhculting the derivation of analytical results. However, a numerical search was made 
for different values of A and some results are shown in Fig. HKb) (for A = 0.6), and (c) (for A = 0.8). An analogy 
between the types of solution for A = 0.5 and solution types for 1/2 < A < 1 can be drawn, having in mind some 
considerations. First, there is no longer a symmetric solution: there is only one solution —the connected solution— 
that exists for Itn —0 and is continuous for all values of Im- Furthermore, instead of appearing or disappearing a 
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FIG. 3: Solutions for A = 1/2. Full curves were used for showing stable solutions, and dotted curves for unstable solutions. 
A number N = 1 + Up + 3ns was assigned to each solution (see caption of Fig. [5] for details). Used parameters {B, aO) are: 
iV = 1, (0.4,1); N = 2, (1.2,1); N = 3, (0.9,1); iV = 4, (0.4, 2.5); = 5, (1.2, 2.5); N = 6, (0.6, 3); N = 8, (1.8, 2.5); = 9, 

(0.5, 7.5). 


couple of asymmetric solutions (as happened for A = 1/2 through a pitchfork bifurcation), now a couple of unconnected 
solutions appear or disappear through saddle-node bifurcations. For —!> 0, if A = 1/2, only the symmetric solution 
exists, if A > 1/2, only the connected solution exists. 

All posible cases can be organized in a simple way similar to what we did for A = 1/2. We define as the number 
of values of where there is a saddle-node bifurcation that makes a couple of unconnected solutions appear or 
disappear (these bifurcations tend to pitchfork bifurcations in the limit A —>■ 1/2) and Ug as the number of regions 
where bistability occurs due to a couple saddle-node bifurcations. 

The possibilities are: only one connected solution (n^ = 0); the connected solution and a couple of unconnected 
solutions between two values of incident field, produced by two saddle-node bifurcations (n^ = 2); and the connected 
solution plus a couple of unconnected solutions starting from an incident field value, produced by a single saddle-node 
bifurcation (n'p = 1). Independently of the value of n^, there may be 0, 1 or 2 regions where bistability emerges due 
to a couple of saddle-node bifurcations on the connected branch (ug = 0, 1 or 2, respectively), or there could be two 
pairs of saddle-node bifurcations on the connected branch and one on an unconnected branch (ug = 3). In order to 
enumerate the cases with a single index, we define = I -I- -I- Sug. No cases have been found with Ug = 3 and 

n'p = 0 . 































































FIG. 4: Examples of different solutions for 1/2 < A < 1. Full (dotted) line was used for showing stable (unstable) solutions. 
The number N = 1 + n'p + 3ns was assigned to each solution, where Us is the amount of couples of saddle-node bifurcations, 
and n'p is the amount of saddle-node bifurcations where an unconnected couple of solutions start or end. Used parameters 
{B, aO, A) were: iV = 1, (0.4,1, 0.7); iV = 2, (1.8, 1, 0.7); N = 3, (1.5, 1, 0.7); iV = 4, (0.4, 3,0.7); iV = 5, (1.4, 2,0.55); N = 6, 
(0.8, 3, 0.6); N = 7, (0.6, 5, 0.7); N = 8, (1.7, 2.2, 0.7); = 9, (1,4, 0.7); iV = 11, (1.8, 4,0.7); N = 12, (1.4, 4, 0.7). 

There is one case in which there is some ambiguity in the above classification: = 1 (unconnected branches) and 

ris = 2 (two regions of bistability). This case corresponds to two qualitatively different solutions: one in which the 
two bistability regions are on the connected solution (this is the case shown in Fig. IH solution 8) and the other in 
which one bistability region is on the connected solution and the other solution is on an unconnected branch. We 
did not find an equivalent case with = 2 and = 2 where one of the bistability regions is on the lower branch of 
solutions. 

Examples of typical solutions are presented in Figured 

There is a last relevant point about solutions with A>l/2: ifnp>la stable unconnected solution with /+ < /_ 
may occur. This is a counterintuitive result, since the pumping in the ‘-I-’ polarization is greater than the one in the 
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’ polarization. An example of this phenomenon is analyzed in the next section. 


IV. NUMERICAL INTEGRATION RESULTS 

In this section, the temporal evolution of homogeneous solutions of Eq. (IT^ is studied. In order to do that, numerical 
simulations applying 4th order Runge-Kutta method are used. 

We found that unstable regions in graphs are due to three possible reasons. Whenever there is a pair of 

saddle-node bifurcations, the middle branch will be unstable, and the other two are expected to be stable, as can be 
seen in, for example. Fig. |3] for V = 4, or Fig. |4] for V = 4. If there is a pitchfork bifurcation (or saddle node that 
generates unconnected solutions), the symmetric branch (or the unconnected middle branch, i.e., the branch closest 
to the connected solution) will be unstable, and it is expected that the other two solutions are stable, as can be seen 
in, for example. Fig. [3] for N = 2 and 3 or Fig. |3]for the same values of N. Finally, there may be areas where we 
expect a particular branch to be stable, but the solution becomes unstable due to a Hopf bifurcation. Examples of 
regions where Hopf bifurcations occur can be seen in Fig. |3]for V = 8 and 9, and Fig. 2] for V = 8, 9, 11 and 12. 

Numerical results show that, if for a given value of there is only one unstable solution, due to a couple of saddle- 
node bifurcations or a pitchfork bifurcation (or saddle-node that generates unconnected solutions), small perturbations 
in the solution make it evolve to one of the two stable solutions, and the chosen solution depends only on the 
disturbance. Also, when there is only one solution and it is unstable due to a Hopf bifurcation, a disturbance takes 
the system to the oscillatory solution. The latter case occurs, for example, with the solution presented in Fig. 2] for 
N = 8 and between 3 and 7.2. Fig. [5] shows the same case, indicating with arrows the possible values of /+ and 
/_ that the unstable solution takes after it is perturbed. It is important to note that the Hopf oscillation occurs in a 
four-dimensional space generated by the real and imaginary parts of A_|_ and A_, and the present graph shows only 
the square modules of those fields. 

Interesting behaviors occur when, for a given value, there are two or more unstable solutions, due to bifurcations 
generated by instabilities of the same or different types. For instance, it can be seen in Fig. [5l for = 20, that 
the unstable unconnected solution evolves towards the stable unconnected solution or to the closest stable connected 
solution (but it does not evolve to the other stable connected solution). Also, in the same figure and for the same 
value of lin, the unstable connected solution (the dotted part of the connected solution, which occurs due to a couple 
of saddle node bifurcations on the connected solution) evolves towards one of the two stable connected solutions (it 
does not go towards the unconnected stable solution). 


V. CONCLUSIONS 

We described the evolution of the electromagnetic field in a ring cavity with flat mirrors, filled with a Kerr-type 
electric and magnetic nonlinear material with positive or negative refraction index. We took into account the degree 
of freedom of the polarization of the incident field. Starting from two pairs of coupled nonlinear Schrodinger equations 
(for the circularly left and right polarized components of the electric and magnetic fields), we found a proportionality 
relation between electric and magnetic field components, so that the description was reduced to just two Lugiato 
Lefever equations. 

Considering that the effective magnetic response can be nonlinear, and taking into account the polarization degree 
of freedom of the incident beam, a rich variety of possible bifurcations in homogeneous stationary solutions was found, 
depending on the values of the constants that describe the system. These solutions were classified and exemplified. 
The number of qualitatively different solutions that were obtained is large, and our classification in terms of the 
number of possible bifurcations simplifies the description. 

In the case of linear polarization (A = 1/2), we found that the solutions /+ and /_ depending on 7^, could present 
between 0 and 2 pitchfork bifurcations and between 0 and 2 pairs of saddle-node bifurcations, yielding 8 possible 
cases. We determined the boundaries between regions of each type of solution in the parameter space given by a© 
and B. The same analysis was performed for circular polarization (A = 1). 

For elliptical polarization (1/2 < A < 1), we found 11 possible cases based on the number of bifurcations. We 
numerically explored regions of parameters that generated each case. 

We also numerically analyzed how unstable homogeneous solutions evolve after a disturbance, and found regions 
where Hopf instabilities take place. 

Finally, for elliptical polarization there are some sets of parameters where, taking an incident field polarized mostly 
in the ‘-I-’ component, it is possible to obtain stable solutions where the ’ polarization component prevails. 

In summary, we presented an analysis of the possible homogeneous stationary solutions of the electromagnetic field 
in a ring cavity containing a material with electric and magnetic nonlinearities, and taking into account the polarization 
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FIG. 5: Example of evolution of unstable solutions (parameters as in Fig. IH = 8). It can be seen that, for the stable 
unconnected solution, 7_ > /+ although the input field has > \Ain-f- Full (dotted) line is for stable (unstable) 

stationary solutions, and gray regions are for oscillations due to a Hopf bifurcation. Arrows show the evolution of a perturbed 
unstable solution for some typical cases. For Ii„ = 2.5, the unstable solution related to a couple of saddle-node bifurcations 
—UCl— evolves into one of the two stable solutions —SCI or SC2—. For /i„ = 9, if the system starts at the connected 
unstable solution —H—, it evolves towards the Hopf oscillation; if the system starts at the unconnected unstable solution 
—UUl— , it can evolve to the lower stable solution —SUl— or to the Hopf bifurcation —H—. For I in = 20, perturbations 
around the unstable point of the connected solution —UC2— make the system evolve to the stable points of the same connected 
solution —SC3 or SC4—; the unconnected unstable branch —UU2— evolves towards the stable unconnected branch —SU2— 
or towards the closest stable point of the connected solution —SC3, it does not evolve towards SC4—. 


degree of freedom of the light. Let us note that even for linearly polarized incident light, elliptically polarized solutions 
can appear, so that considering only a scalar equation may cause the removal of a possible solution. Despite the large 
number of parameters of the original description, a classification of the possible solutions for a given polarization could 
be obtained taking into account only two parameters: a (related to the sign of the electric and magnetic nonlinear 
susceptibilities) times 0 (the scaled detuning) and B (related to the electric and magnetic nonlinear constants Be 
and Bm)- Results presented here can be useful when ring cavities are used in optical switching or information storage. 
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